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Topological quantum computation has been extensively studied due to its robustness against 
decoherence. A conventional way to realize it is by adiabatic operations—it requires relatively long 
time to accomplish so that the speed of quantum computation slows down. In this work, we present a 
method to realize topological quantum computation by periodic driving. Compared to the adiabatic 
evolution, the total operation time can be regulated arbitrarily by the amplitude and frequency of 
the periodic driving. For the sinusoidal driving, we give an expression for the total operation time 
in the high-frequency limit. For the square wave driving, we derive an exact analytical expression 
for the evolution operator without any approximations, and show that the amplitude and frequency 
of driving field depend on its period and total operation time. This could provide a new direction 
in regulations of the operation time in topological quantum computation. 


I. INTRODUCTION 


Decoherence is an enemy of quantum computation, 
which is the loss of coherence due to the presence of en¬ 
vironments. As a promising avenue to deal with the de¬ 
coherence, topological quantum computations [l|-[8| em¬ 
ploy two-dimensional quasiparticles called anyons, whose 
world lines cross over one another to form braids in a 
three-dimensional spacetime. Information encoded in the 
anyons is robust against local perturbations and quan¬ 
tum operations can be performed by braiding the non- 
Abelian anyons Hi- The simplest example of the non- 


Abelian anyons is the Major ana fermions which are 
dieted to exist in fractional quantum Hall systems 


Dre- 


topological insulators 


tionalqn 

00 , 


solid state systems 


10 |, 

13 , 

14- 


and semiconductor-superconductor hybrid systems ^ _ 

The signatures of Major ana fermions have also been 
observed in experiments more recently [Siil , which 
gives rise to an opportunity to encode a qubit by Ma- 
jorana fermions in these materials. 


A quantum task is often accomplished by a sequence 
of quantum operations rather than single quantum oper¬ 
ation [ 2214 ^ . The total operation time increases linearly 
with the increasing of the number of operations. Consid¬ 
ering the limit of the coherence time of quantum systems, 
long operation time is not favorable, even if the quantum 
topological computation is robust against perturbations. 
On the other hand, the time-periodic driving systems 
have been extensively studied in the past few years. Es¬ 
pecially, several work HiEl have shown that the topo¬ 
logical properties can be changed in topologically triv¬ 
ial system by time-periodic driving (e.g., the existence 
of Floquet topological insulators or Floquet Majorana 
fermions). Recently, the Floquet Majorana fermions is 
realized by periodic driving fields in the system of cou¬ 


pled quantum dots proximity to a s-wave superconductor 
0 . More recently, it has been proposed to achieve the 
direct coupling between the top olog ical and conventional 
qubits by periodic driving fields Q . In this paper we ex¬ 
plore the possibility to regulate the total operation time 
by periodic driving For concreteness, the physical model 
of interest is the quantum dots coupled to the Majorana 
modes in a topological superconductor. Of course, this 
method can also be extended to the other quantum sys¬ 
tems. 

The paper is organized as follow. In Sect. nn we briefly 
introduce the topological quantum computation by adi¬ 
abatic evolution. In Sect. uni we first recall the Floquet 
theory, then the periodic driving fields in the form of si¬ 
nusoidal, square wave, and (5-function kick are applied 
separately to modulate the total operation time for re¬ 
alizing the quantum operations. Finally we extend this 
method to other hybrid quantum systems in Sect. IIVI 
The discussion and conclusion are given in Sect. [Vl 


II. QUANTUM COMPUTATION BY 
ADIABATICAL EVOLUTION 


Recently, the adiabatic evolution has widely applied to 
the preparation and manipulation of Majorana fermions 
1^ . 1^, 0. In particular, it has been shown that topo¬ 
logical quantum information processing becomes possible 
in the one-dimensional network by adiabatically con¬ 
trolling the locally tunable gates which affect the chemi¬ 
cal potential over a finite length of the wire. In following 
we describe the main idea of adiabatic evolution. That 
is, design a Hamiltonian Hi whose ground state is the 
target state |4 ^t) while the ground state |4^o) of Hamil¬ 
tonian Ho is easily to prepared. Assume that there exists 
a quantum system satisfying the following Hamiltonian 
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where /(t) is a slowly varying function of evolution time 
t with /(O) = 0 and /(I) = 1. According to the adia¬ 
batic theorem, the quantum system evolves adiabatically 
from the initial (ground) state |^o) to the target (ground) 
state I^t) at time t = T. 

In the present work, the physical model of interest 
consists of a quantum dot coupled to a semiconducting 
nanowire, as shown in Fig. HJa). In a magnetic field, by 
proper spin orbit interaction and proximity coupling to 
a superconductor, the nanowire can exist the Major ana 
bound states in the topological phase [H EE El m. 
Then the effective Hamiltonian (in the low-energy limit) 
for th e q uantum dot coupling to the Majorana mode 
reads [ij 

H = s{t)d^d -h {v'^d^ — vd)^i, (2) 

where a (a’*’) is the annihilation (creation) operator for 
the electron in quantum dot and the on-site energy e{t) 
for the quantum dot can be controlled by the gate volt¬ 
age Vg. V denotes the tunnel coupling between the quan¬ 
tum dot and the Majorana mode 71 . Without loss of 
generality we assume v is an real number and take all 
physical parameters in units of v. Since the Majorana 
mode is Hermitian ( 7 J = 7 ^ and 7 ^^ = 1), we cannot 
use the number operator 7 j 7 i to count the occupation 
of the Majorana mode. Whereas, two Majorana modes 
can be combined to generate one ordinary fermion, e.g., 

and 72 = i{b^ — b). One can adopt the num¬ 
ber operator b^b of the ordinary fermion to count the 
Majorana modes. 



FIG. 1: The setup for realizing the operation Pi. 


Since the total parity of the electron in quantum dot 
and the ordinary fermion formed by Majorana modes is 
conserved, the Hamiltonian is block diagonal in the basis 
spanned by {|0)F|0)i:), \1)f\'^)d, |l)F|0)i:), |0)F|l)i:>}, 


/ 0 V 0 0 \ 

V e{t) 0 0 I 

0 0 0 ^ ’ 

V 0 0 V e{t) ) 


(3) 


where the state \m)F\n)D {m^n = 0 , 1 ) represents m or¬ 
dinary fermion formed by Majorana mode and n electron 
in the quantum dot. In Ref. [^, it suggests that by adi¬ 
abatically changing the values of ^ from — oc to +oc, it 
can realize the operation Pi which denotes the inversion 


of the occupation in ordinary fermion combined by the 
Majorana mode (i.e.. Pi = 71 ), 

Pi(sin6>|0) F + cos 6 >|l) f) = sin 6 >|l) f + cos0\0)f‘ (4) 

Fig. [ 2 ] shows the different dynamics behaviors for dis¬ 
tinct changing rate It can be observed in Fig. [Sfa) 
that the operation Pi cannot be achieved perfectly since 
the changing of ^ does not satisfy the adiabatic condition 
very well (It cannot satisfy H(t) <C 1 all the time). Thus 
the changing rate ^ should be small in order to meet the 
adiabatic condition, along with the increasing of opera¬ 
tion time (cf. Fig. [2fc)-(d)). To implement the single 
qubit rotation or non-Abelian operation, one shall suc¬ 
cessively execute the operation Pi. Therefore the total 
operation time increases with the increasing of the num¬ 
ber of operation Pi. In addition, it needs to point out 
that this Majorana based qubits may be susceptible to 
decoherence due to the electron tunnel coupling process 
liiiil- As a consequence, the adiabatic evolution is at a 
disadvantage in minimizing the influence of decoherence 
as far as possible. Recently, it is overcome by short¬ 
cuts to adiabaticity for the non-Abelian braiding with 
Y-junction structure [52|. In following it demonstrates 
that the situation can be changed by periodic driving. 



FIG. 2: The time evolution of the system to obtain the oper¬ 
ation Pi, where the fidelity is defined as |('0(t)|'0(T))|^. The 

expression of V{t) is defined as, V{t) = I I where 

\Ei{t)){i = 1,2) is the instantaneous eigenstate with corre¬ 
sponding eigenvalue Ei{t) in the even (or odd) subspace (The 
adiabatic condition can be written as V{t) <C 1). The initial 
state is |'0(O)) = sin^|0)F + cos^|1)f and the target state is 
|'0(T)) = sin^jl)F + cos^|0)f, 0 — We have set the on-site 
energy of the quantum dot increases with time linearly, e.g., 
£(t) = The operation time is 2T and the final value of 
e{t) is 50 during the time evolution. All parameters are in 
units of the tunnel coupling v. (a) T = 10. (b) T = 30. (c) 
T = 50. (d) T = 80. The larger T means the smaller of the 
changing rate A One can find that the operation time for 
perfectly achieving the operation Pi is about 2T = 100. 
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III. QUANTUM GATES WITH PERIODIC 
DRIVING 

A. Floquet theory 

Let us first recall the Floquet theory briefly [53| . Pro¬ 
vided that the system Hamiltonian has a time-periodic 
driving field, H{t) = H{T + t), where T is the period 
and the driving frequency reads uj = ^. The Floquet 
theory asserts that the solutions of Schrodinger equa¬ 
tion have the form |^n(^)) = is quasi¬ 

energy and the Floquet state |0n(^)) has the property 
\4^n{t)) = |0n(F + t)). They are satisfied the following 
eigenvalue equation [h = 1) 

=en\(l)„{t)), (5) 

where Heff = is defined as the Floquet Hamil¬ 

tonian. To solve this equation it is very instructive to 
introduce an extend Hilbert space [HJ of time-periodic 
functions with the inner product (('I*)) = ^ dt{-\-). 

With regard to the periodic driving system, it is nec¬ 
essary to make definite on the time-scales during the 
evolution. For the case of Floquet state |0n(^))5 since 
it has the same period with the driving field, it affects 
the system dynamics on short time-scale (in the high- 
frequency limit). What really affects the long time-scale 
of the system dynamics is the gap of the quasi-energies. 
Therefore, it is crucial to determine the evolution time 
through modulating the structure of quasi-energies in the 
periodic driving system. 


B. Sinusoidal driving 

We first consider the periodic modulation of the on¬ 
site energy for quantum dot with the sinusoidal form 
e{t) = 5ocos(co’t), which can be created by a waveform 
generator. In order to obtain an approximate expres¬ 
sion for the quasi-energy, we solve the time-dependent 
Schrodinger equation by standard perturbation theory 
[lEii, where the tunneling Hamiltonian is regarded as 
the perturbation. Due to the total parity conservation it 
is convenient to study in the even parity (or odd parity) 
subspace. Then the Hamiltonian reduces 2x2 matrix. 
Since iLe(t) = e{t)a^a is diagonal, when substituting into 
Eq.([5]), the eigenstates of can be readily 

given by 


|Ai(t)) = (e^"^^0)^ 
|A2(t)) = (0,e' 


iX2t—i— sin ujt\T 


(6) 


where Xi (z = 1,2) is the corresponding eigenvalue (i.e., 
quasi-energy). On the other hand, due to the period 
of Floquet states, one then can find that the zeroth or¬ 
der approximation of both quasi-energies are zero (mod¬ 
ulo uj). Thus the time-dependent eigenstates can be 


approximately viewed as time-independent eigenstates 
|Ai(t)) = (1,0)^ and |A2(t)) ^ (0^1)^ in the high- 
frequency limit (cj ^ 1). The first order approximation 
of quasi-energies can be obtained via diagonalizing the 
perturbing matrix [56| 


?). ( 7 ) 

where the matrix element Q = ^ af¬ 

ter some straightforward calculations. Consequently, the 
quasi-energies are calculated as ei^2 = i|Q| (in the “first 
Brillouin zone”) and the corresponding eigenstates be¬ 
come |ei,2) = :^(|Ai(t)) ± |A2(t))). The gap of quasi¬ 
energies are then given hy A = 2\Q\. In the light of the 
identity 


sin ujt 




(8) 


where J'n is the n-order Bessel function, we can finally ob¬ 
tain the analytical expression for the quasi-energies gap 
A = 2\vM^)\. 


<1 1 


o'-^^^^— - 

0 20 40 ^ 60 80 100 



FIG. 3: The gap of quasi-energies versus (a) the amplitude eo 
when ct; = 10, (b) the frequency cu when eq — 50. The gap of 
quasi-energies approaches 2 when the driving frequency tends 
to 60 in panel (b). After that the gap increases slowly with 
the increasing of driving frequency. The dynamics evolution 
of periodic driving system with (c) uj = 10, (d) uj = 20.8, (e) 
Ct; = 50. All parameters are in units of the tunnel coupling v. 

Fig. [31(a)- (b) demonstrate the relation between the 
driving field and the quasi-energies gap, while Fig. |3](c)- 
(e) show that the evolution time reduces with the increas¬ 
ing of the quasi-energies gap. Therefore, we can choose 
special evolution time for the periodic driving system by 
appropriately selecting the frequency and amplitude of 
driving field. An inspection of Fig. |3Ka) also shows 
that the operation time varies with the decreasing of the 
amplitude of driving field when we fix a high frequency. 
Interestingly, there exists a special case that the quasi¬ 
energies of periodic driving system vanishes (namely the 
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two quasi-energies approach degeneracy) by choosing the 
amplitude Eq and the driving frequency uj properly, which 
is known as coherent destruction of tunneling [^, . As 

a consequence the state is localization so that it is invalid 
to achieve the operation Pi, as shown in Fig. (H^d). 

Since the initial state can be approximately given by 
I'0(0)) ^ ^(|ei) — 02)), the time evolution of periodic 
driving system approximately reads 

m)) ^ ^(e--^*|ei)-e--^‘|e2)) 

~ cos[Jo{ — )vt]\X 2 ) -isin[J'o(—)i'i]|Ai).(9) 

UJ UJ 

Hence the total operation time for realizing the operation 
Pi approximately equals 




( 10 ) 



FIG. 4: The total operation time T versus the amplitude eo 
and the frequency uj of the driving field. 

Fig. 0] depicts the relation between the total operation 
time and the amplitude as well as the frequency of the 
driving field. It suggests that one shall avoid the param¬ 
eter regions with the coherent destruction of tunneling, 
since it takes long operation time to realize the operation 
Pi. Apart from this regions, the total operation time can 
be regulated within proper range. In addition, from the 
Eq. m, one can find readily that the expression of fi¬ 
delity is 

F = sin2 [Joi-)vt] = hi - cos [2M^)vt]). (11) 

UJ Z UJ 

Fig. [5]plots the relation between the gap of quasi-energies 
and the coefficients of fidelity in the exact and pertur¬ 
bation regime, respectively. It demonstrates that the 
perturbation results work extremely well in the high- 
frequency limit. 

In order to check the validity of the perturbation the¬ 
ory, we plot the dynamics of the system with different 
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FIG. 5: The coefficients (a) ao, (b) ai, (c) bi, (d) u' of fidelity 
versus the gap of quasi-energies, where the expression of fi¬ 
delity for curve-fitting is P' = ao + cn cos uj't + bi sin cjT. The 
lines represent the analytical solutions given by Eq. cu with 
perturbation theory while the circles, squares, and stars rep¬ 
resent exact results obtained by curve-fitting. Note that the 
curve-fitting has high degree of precision for the exact results 
since the values of R-square and Adjusted R-square approach 
unit (> 99.37%) in MATLAB. 


driving frequencies. The results are given in Fig. [6l We 
observe that the dynamics is well in agreement with the 
results by perturbation theory when u; > 10 (in units of 
u), while it deviates seriously from the perturbation re¬ 
sults when cj < 10 (in units of see the pink dot-dash 
line in Fig. [6j). As a result, one can employ the perturba¬ 
tion theory safely when the frequency of the driving field 
is at least an order of magnitude larger than the tunnel 
coupling. 



FIG. 6: The evolution of the fidelity with different frequencies 
of driving field. We have set the quasi-energies gap A = 
2\vJq{^)\ — 1.0236, thus the period of system dynamics is 
approximated to T ~ ^ = 6.1383 in the high-frequency limit, 
which is confirmed by the green solid line and red dash line. 

As mentioned, the perturbation theory is not valid in 
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the low-frequency limit, this gives rise to a question how 
the system behaviors in this limit. Now we go to explore 
this issue. When the driving frequency is small such that 
the adiabatic condition approximately holds (since the 
on-site energy changes slowly), we expected that the sys¬ 
tem dynamics, e.g., the fidelity, at the long time-scale 
would be periodic with period T = —. As expectation, 
we find from the dash line in Fig. [71^a)-(c) that this is 
exact the case. Besides, one can observe that the high 
fidelity lasts a long time within a period when the driv¬ 
ing frequency is small, see Fig. [7fa)-(d). The fidelity 
changes fast when the frequency of driving field is large 
(see Fig. [7fa), where the yellow lines and blue lines alter 
frequently). It also affects the dynamics when the ampli¬ 
tude of driving field is large, which is shown in Fig. [Tl^b). 
Fig. [71^c) illustrates how the offset energy affects the fi¬ 
delity. Interestingly, the high fidelity lasts longer time 
(see the yellow region) when the offset on-site energy is 
larger. 



FIG. 7: The fidelity as a function of the evolution time and 
the distinct parameters of driving field. The expression of 
driving field is e{t) — sq cos cut —where ^ is the offset energy 
of driving field, (a) eo = 40,.f = 20. (b) uj = 0.05,^ = 0.5£o. 
(c) eo = 40,ct; = 0.05, = (5£o. (d) eo — 40,uj = 0.04, = 20. 

The dash lines are ploted by the function t — where n 

is integer. 


C. Square wave driving 

It is believed that the periodic square wave driving 
fields are easily achieved in practice. In fact these driv¬ 
ing fields have been studied extensively in time-periodic 
driving system. In particular, it has been shown in ex¬ 
periment [Hil that the Stiickelberg interference in a su¬ 
perconducting qubit is driven by the square wave form 
which we use in following. We first give the exact ana¬ 
lytical expressions for system evolution operator without 
any approximations. The square wave driving for the 


on-site energy is expressed as 


f 5i, nT < t < ti + nT, 
[ ^2, ti + fiT < t < (n -j- 1)T, 


( 12 ) 


where n = and ti G [0,T]. Thus the system 

evolution operator U of one period can be written as 
U{T^0) = with t 2 = T — ti. After some 

lengthy algebra, one can obtain the expression of the evo¬ 
lution operator 


U{T,0) = 




X 1 X 2 


A B^C 
B-C A 


(13) 


where 
A = 


A 2 • X\t\ , X2^2 / X\t\ , . X\t\ . 

—Av sm-sm-h [xi cos-h isi sm-) 

2 2 ^ 2 2 ^ 

/ X2t2 , . . ^ 2^2 X 

\X2 COS —--h 182 sm 


o . X2t2 . Xiti Xih . X2t2. 

B = — 2 iv(x 2 COS -sm-h Xi cos-sm-), 

^ 2 2 2 2 ^ 

C = 2u(5i- 52)sm^—sm^—, 


D = eiti^S2t2, xi = X2 = Jel + Av‘^. (14) 


At first we design the driving time ti (^ 2 ) of the on-site 
energy Si ( 62 ) to satisfy XiB = tt (^ 2^2 = tt), that is. 


h = 


TT 

a/4u 2 -h sj ’ 


t 2 


TT 


tl + t 2 


T. (15) 


Consequently, the period of the square wave driving is 
confirmed. According to Eq. (USD, the evolution operator 
can be further simplified. 


U{T, 0) = — -e-i^ ( , (16) 

X1X2 \X4, X3 J ' 


where xs = Av^ -h 51^2 and X 4 = 2v{si — ^ 2 )- After N 
evolution periods, the final evolution operation becomes 


u{r^o) 


y i{ri — r^) 


(17) 


where ri = x^—ix^ = \r\e r 2 = x^A-ix^ = |r|e*^, |r| = 
^^3 + X 4 , and tan 6 > = From the expression in Eq. 

it clearly requires -h = 0 in order to realize 
the operation Pi perfectly (up to a global phase factor). 
By making the vectors and produce a 7 r-phase 
difference, that is, NO — {—NO) = tt, one can readily 
obtain the number of evolution periods 


N = 


TT 


2 arctan 


2i;(£2-£i) ‘ 
Av‘^-\-S\S2 


(18) 


Note that according to Eq. ([T 8 |) the number of evolu¬ 
tion periods N is not integer generally. Nevertheless it 
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does not affect the main results because we can just take 
an integer nearest to N, as the fidelity increases slowly 
when it approaches 1. In turn, if one designates the pe¬ 
riod T and the number of evolution periods N in the 
periodic square wave driving system, the values of on¬ 
site energy ei and 82 can be determined by Eq. m and 
Eq. ([T 8 |) as well. It demonstrates the system evolution 
with distinct values of ei and 82 in Eig. [Hl(a)-(b), as well 
as the special period and the number of evolution periods 
in Eig. [ 8 ](c)-(d). As expected, it can also realize the op¬ 
eration Pi by square wave driving and we can modulate 
the period and the total operation time in this case. 


In absence of (5-function kick, the expression of fidelity 
for realizing the operation Pi becomes sin where 
the maximum of fidelity is | ^ P • One easily observes that 
it cannot obtain the operation Pi when ^ 0 , i.e., the 

on-site energy 81 ^ 1. However the situation changes in 
the presence of (5-function kick and it can realize the op¬ 
eration Pi regardless of the large value of on-site energy 
81 (the value of 81 only determines the driving period). 
Especially, when we take away the (5-function kick if the 
operation Pi has been completed, the system is still sta¬ 
tionary. 






FIG. 8: The system dynamics of distinct forms for square 
wave, (a) si = 40, 52 = 60. (b) si = 10, 52 = 90. (c) 
T = 0.3, N = 60. (d) T = 0.2, N = 30. The other two 
parameters are calculated by Ea. (p31) and Eg. ITSl) . 


D. (5-function kick 

When 52 ^ + 00 , one can readily find in Eq. ([T 5|) that 
t 2 ^ 0. Then the square wave driving field reduces to 
periodic (5-function kick, i.e.. 


N 


z{t) = 5i -h 52 ^ (5(t - nT), 


(19) 


n=l 


where T is the driving period and can be calculated ap¬ 
proximately as T ^ Therefore the total oper- 

ation time for realizing the operation Pi is approxima- 
tively 


r ^ 


TT 


2 ^JAv‘^ + 8 \ arctan 


( 20 ) 


Note that the dynamics behavior is quite different from 
the absence of ^-function kick, i.e., under a static driving 
field. In the static case, the evolution operator U reads 


l/(t, 0 ) = e 


i£\t 


= \f^ 


COS ^ -h sin ^ 

^ • 2 ?; ^ Tf ^ 

sm^ 


sin f 


sin ^ 


IV. APPLICATION TO OTHER SYSTEMS 

The periodic driving method can be applied to the 
other structure of hybrid quantum dot-topological sys¬ 
tem. Here we apply it to a system described by the fol¬ 
lowing Hamiltonian [ 6 l|, as illustrated in Eig. 

R = 5(t)(d|d>^ -h djaj,) + Ed|d>^d|d^ -h — via^)^i 

+(t>2a| - V2ax)^2- (22) 

8{t) is the on-site energy of the quantum dot. Vi{i = 1 , 2 ) 
denotes the tunnel coupling between the quantum dot 
and the Majorana mode 7 ^. In particular, the spin-up 
(labeled as t) and spin-down (labeled as |) electrons can 
only tunnel into the Majorana mode 71 and 72 , respec¬ 
tively. V represents the energy contributed by double 
occupation on the quantum dot. In the situation of large 
V, the quantum dot can only hold single electron. 

Since the total parity (the electrons in quantum 
dot and the ordinary fermions formed by Majorana 
modes) of the hybrid system is conserved, we can 
restrict ourself in the even-parity subspace spanned 
by {|0 )fi|0)f2|0)d, |1)fi|1)f2|0)d, |0)Fi|l)F2|lt)D> 

|0)Fi|l)F2|li)D, |l)Fi|0)F2|lt)F>, |l)Fi|0)F2|li)D,}, 

where the subscript Fi(i = 1 , 2 ) represents the ordinary 
fermions formed by Majorana modes. The matrix form 
of Hamiltonian then can be written as 


H = 



0 

0 

V2 

Vi 

0 \ 

0 

0 

Vi 

0 

0 

V2 

0 

^1 

8(t) 

0 

0 

0 

^2 

0 

0 

5(t) 

0 

0 


0 

0 

0 

5(t) 

0 

0 

^2 

0 

0 

0 

8 {t) ) 


(23) 


-h 4v^. 


This quantum dot-Majorana system [60|, can be 
used to prepare entanglement between spin and topologi¬ 
cal qubits or quantum information transfer between spin 
and the topological qubits (even for the quantum logic 
gates) by the adiabatic evolution. We take the prepa¬ 
ration of entanglement (denoting as the operation P 2 ) 
between the electron spin and Majorana modes as an ex¬ 
ample to exemplify how to manipulate the operation time 
(21) by periodic square wave driving given in Eq. ([12]). The 
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FIG. 9: The setup for realizing entanglement between the 
convention qubit and topological qubit. 


operation P2 reads, 

P 2 (sin6>|0)i7j0)i?2 + cos6>|l)i7jl)i72)|0)i:> 

= ^(sin6>|l)i7jO)F2 + cos6>|0)Fjl)F2)|lt)D 

Vo 

H-(sin6>|0)Fi |1 )f2 + cos6>|1)fi |0)F2)|lt)^-(24) 

V 

where v = As the Hamiltonian is a 6 x 6 ma¬ 

trix, the analytical expression of the evolution operator 
U{T^ 0 ) = involved. Here we only give 

the equations that determine the period of the driving 
field and the total number of evolution periods, i.e.. 


T — ti + ^2, 

^ ( 25 ) 

2 arctan -4- 


where p 


yGf+ 4 ( 7 f+ 7 |) 


, h = 


yj|+ 4 ( 7 f+ 7 |) ’ 


Xo = 


S 1 S 2 + ^vl + Avl, and x'^ = 2(^2 - Fig. 

[To] plots the fidelity of realizing operation P2 as a func¬ 
tion of the evolution time by the adiabatic evolution and 
the periodic square wave driving, respectively. Again, 
we find that the operation time for adiabatic evolution 
requires relatively long time since it must satisfy the adi¬ 
abatic condition while the operation time and the period 
of the square wave driving can be regulated. 



FIG. 10: (a) Realizing the operation P 2 by the adiabatic 

evolution. The on-site energy of the quantum dot increases 
with time slowly, e{t) — 1.25t. (b) Realizing the operation P 2 
by the periodic square-wave driving, si — 30, 82 — 50. 0 — 
r’l = r ’2 = All parameters are in units of v. 


V. DISCUSSION AND CONCLUSION 

In the last section, we have studied how to implement 
the operation Pi by periodically driving the on-site en¬ 
ergy of the quantum dot. For a single operation Pi, it is 
far from sufficient to permit quantum computation. We 
next briefly discuss how to realize an arbitrary rotation 
for a qubit by successively executing the operation Pi 
twice. 

As shown in Fig. [Hb), the system Hamiltonian of in¬ 
terest reads 

H = e{t)o^a-\-{\vi\e~'^^^a^ — 

+(|v 2 |e“*'^^a^ - |v2|e*®"a)72, (26) 

where we have introduced a phase 0i {O2) into the tunnel 
coupling vi (V2) in the Hamiltonian. Defining the opera¬ 
tor c = 2: = ^2 = 

the Hamiltonian (| 26 |) becomes, 

H = £{t)8c + z{c^ - c){ziji + 2 : 272 ), (27) 

where the phase difference 0i — O2 equals 2 niT {n is integer 
and the phase can be modulated by the magnetic flux 0). 
The form of Eq. (| 27 |) is the same as Eq. ^ if we redefine a 
new Majorana mode 712 = ^171 + ^272, where the tunnel 
coupling is denoted by Clearly, the operation Pi = 712 
in this notation. Consider the two level system spanned 
by {|0)i7, |l)i7}, then we can express the Majorana opera¬ 
tors in terms of Pauli matrices cr^, cr;^}, i.e., 71 =0-3,, 
72 = cFy, 7172 = By successively executing the op¬ 
eration Pi twice with different relative tunnel coupling 
strengthes between modes zi and 2(2, the total operation 
becomes P = 7127I2 = (^i 7 i + ^ 272 )( 47 i + ^72) = 
{ziz[ Z2Z2) i{ziz2 — ZiZ2)cFzi which is exactly an ar¬ 
bitrary rotation around the 2(-axis. 

Due to the conservation of the total parity, a qubit 
shall be encoded by four Majorana modes The 

Majorana-based qubit can be realized by the general¬ 
ization model in Eig. [IJb), which describes a system 
consisting of three quantum dots coupling to four Ma¬ 
jorana modes ( 71 , 72, 73, 74) in the topological super¬ 
conductor with comb structure. In the even-parity sub¬ 
space spanned by {|0)i7j0)i72 5 \'^)Fi\'^)F2}^ fbe operation 
P = 7i27i 2 is in fact the rotation around the z-axis, and 
the operation P = 723723 is the rotation around the pr¬ 
axis (723 = Z242 + ^373)5 where the ordinary fermion Pi 
(P2) is formed by the Majorana modes 71 and 72 (73 and 
74 )- 

Generally speaking, the tunnel coupling between the 
quantum dot and the Majorana mode depends on both 
the differences among the on-site energies and the tun¬ 
nel barriers. By making use of the periodic driving on 
the on-site energy of the quantum dot, the tunnel cou¬ 
pling would change consequently. Reminding that we can 
employ additional electrostatic gates to manipulate the 
tunnel barriers, the tunnel coupling can keep a constant 
in practise, even the on-site energies change. Indeed, the 
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possibility of controlling the tunnel coupling in semicon¬ 
ductor nanowire has been experimentally shown recently 
[63|. So, we believe that in the near future it is also possi¬ 
ble to manipulate the gates such that the tunnelling rate 
remains unchanged in our case, especially with the peri¬ 
odic square pulses (since it has only two distinct values 
of on-site energy). 

In conclusion, we have proposed a method to regu¬ 
late the total operation time of quantum computation, 
which can be achieved by periodic driving. By solving 
the time-dependent Schrodinger equation with perturba¬ 
tion expansion, we have given an expression of the quasi¬ 
energies and elucidated the relationship between the total 
operation time and quasi-energies in the high-frequency 
limit. As a result, the operation time can be manipu¬ 
lated by designing the amplitude and frequency of the 
driving field. For the case of low-frequency limit, due to 
the invalidity of the perturbation theory, we study the 
dynamical behaviors by numerical simulations. We find 
the results approach those given by adiabatic evolution. 
Different from the adiabatic evolution, the system in the 
low-frequency limit manifests more intricate behaviors 
and the operation time can also be regulated by the driv¬ 


ing field. In particular, the total time that the high fi¬ 
delity lasts are closely related to the frequency and offset 
energy of driving field. For the case of periodic square 
wave driving, we have derived an analytical expression 
for the evolution operator without any approximations. 
By this expression, we can calculate the amplitude of 
square wave driving fields with fixed operation time and 
period of driving fields. We have also discussed the real¬ 
ization of quantum operations by the (5-kick, which can 
be treated as a deformed square wave driving. The pe¬ 
riodic driving can also be applied to the other quantum 
system—it opens up a new avenue in manipulations of 
operation time in topological quantum computations. 
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